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A generic quantum channel can be represented in terms of a unitary interaction between the information- 
carrying system and a noisy environment. Here, the minimal number of quantum Gaussian environmental 
modes required to provide a unitary dilation of a multi-mode bosonic Gaussian channel is analyzed both for 
mixed and pure environment corresponding to the Stinespring representation. In particular, for the case of pure 
environment we compute this quantity and present an explicit unitary dilation for arbitrary bosonic Gaussian 
channel. These results considerably simplify the characterization of these continuous-variable maps and can be 
applied to address some open issues concerning the transmission of information encoded in bosonic systems. 
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I. INTRODUCTION 

Bosonic Gaussian channels (BGCs) are an important spe- 
cial class of transformations that act on a collection of bosonic 
modes preserving the Gaussian character of any Gaussian in- 
put quantum state flJJ] . The set of BGCs is singled out from 
a theoretical perspective yj], but most significantly also from 
the perspective of practical implementations, since it emerges 
naturally as the fundamental noise model in several experi- 
mental contexts. In the vast majority of physical implemen- 
tations of quantum transmission lines quantum information 
is almost invariably sent using photons — through optical 
fibers yfl, in free space j^D, or via superconducting transmis- 
sion lines J5t]) — physical situations for which BGCs provide 
extraordinarily good models. What is more, BGCs play a ma- 
jor role in characterizing the open quantum system dynamics 
of various setups which use collective degrees of freedom to 
store and manipulate quantum information, including systems 
from cavity QED, nano-mechanical harmonic oscillators [6], 
or clouds of cold atomic gases |7[]. 

Unsurprisingly, therefore, the study of Gaussian or quasi- 
free quantum channels has a long tradition JH 0, d] . Intense 
recent research has mostly been focusing on properties of 
BGCs with respect to their ability to preserve and transmit 
quantum information (for a review see, e.g., Ref. [9] and ref- 
erences therein). Recent contributions include the computa- 
tion of the quantum capacity ifloll of a large subset of single 
mode BGCs 111 111 , a characterization in terms of the notion of 
degradability — introduced in Ref. [ 12] — that allows one to 
identify a zero-quantum capacity subset of BGCs, and a nec- 
essary and sufficient condition for BGCs being entanglement- 
breaking 111 311 . A general unitary dilation theorem for BGCs 
was proven in Ref. 01: It shows that each BGC $ acting on a 
system A formed by n input bosonic modes admits a unitary 
dilation in terms of a bosonic environment E composed of 
I ^ 2n modes, the initial state pe of which is Gaussian, with 
a Gaussian unitary coupling Ua.e corresponding to a Hamil- 
tonian that is quadratic in the canonical coordinates, 

= Tv E [U A , E (p ® p E )U\j ■ (1) 
Here, p is the input quantum state of the system A and Tr# de- 



notes the partial trace over the degrees of freedom associated 
with E. 

The fact that the number of environmental modes i enter- 
ing in the unitary dilation can be bounded from above by 2n 
may be viewed as the continuous-variable counterpart of the 
upper bound on the minimal dilation set by the Stinespring 
theorem lfl6ll for finite dimensional quantum channels: It in- 
dicates that any quantum channel can be described by using 
an environment which is no more than twice the size of the 
input system. An important open question is the characteriza- 
tion of the minimal value of £^ that is needed to represent 
a given BGC, specifically the minimal value £ p *e in case of a 
pure unitary dilation. Similarly to the quantum capacity, this 
quantity may be used to induce a partial ordering in the set of 
BGCs since, as a general rule of thumb, one expects that the 
larger it is, the noisier and the less efficient in preserving the 
initial state will be the associated channel. Furthermore, an 
exact estimation of such number will allow one to consider- 
ably simplify the degradability analysis of BGCs by minimiz- 
ing the number of degrees of freedom of the corresponding 
complementary channel. 

The main result of this work is to explicitly identify this 
minimal value £ p *e — so the minimum number of environ- 
mental modes initially in a pure Gaussian state pg in a unitary 
dilation ([TJ (related to the Stinespring dilation lfl7ID — and to 
construct the corresponding dilation. To simplify terminol- 
ogy, in this case we speak of Eq. (Q~|i as of the Stinespring rep- 
resentation. This is accomplished by first determining a lower 
bound for £^ e in terms of the minimum number g£v of ancil- 
lary modes which are needed to construct a Gaussian purifica- 
tion of a (generalized) Choi-Jamiolkowski (CJ) Gaussian state 
of $. This is motivated by the fact that any Gaussian Stine- 
spring representation (Q~|i naturally induces a Gaussian purifi- 
cation of the CJ states of the channel. Then, we show that this 
lower bound can be exactly achieved by explicitly construct- 
ing a Gaussian Stinespring dilation with q^ n modes. In the 
second part of the paper we finally address the case of unitary 
dilations (Q]) in which the environment state ps is not neces- 
sarily pure, and provide an estimation for the minimal t which 
improves the one presented in Ref. Il4tl . 
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The paper is organized as follows: In Sec.|II]we recall some 
basic definitions and set the notation. The notion of a gener- 
alized CJ state of a BGC and the lower bound for £oure 
are presented in Sec. [Til] Then, in Sec. |IV] we present an 
explicit recipe to construct such minimal dilation. The case 
of dilations involving not necessarily pure environments is fi- 
nally addressed in Sec.|V] while conclusions are presented in 
Sec.[VT] This work includes also some technical appendixes. 



II. BOSONIC GAUSSIAN CHANNELS 

Consider a system A composed of n bosonic quantum me- 
chanical modes described by the canonical coordinates R := 
(Qi, ■ ■ ■ , Q n ] P\i ' " ' j Pn) and by the Weyl (or displacement) 
operators 



V(z) 



(2) 



with z := (xi, ■ ■ ■ ,x n ,yi, ■ ■ ■ ,y n ) T £ R 2 being a column 
vector To simplify notation, we choose units in which 
h = 1. A BGC $ acting on A is completely determined by 
assigning a real vector v G R 2 ™ and two 2n x 2n real ma- 
trices Y, X £ GZ(2n,R) satisfying the complete positivity 
condition 



Y 3? zE with £ := a 2n ~ X T a 2n X . 



(3) 



where a% n is the matrix defining the symplectic form capturing 
the canonical commutation relations of n modes, i.e., 



C2r. 






In " 


-In 






(4) 



with l n indicating the n x n identity matrix. More precisely, 
the map $ is defined as the linear mapping which, for all z 
complex, induces the following transformation 



<t>{p\z) ^ 0($(/5);z) 
where 



±1 ~ v \ —\z T Yz J riv T z 



<j>(p;z) i=Tt[pV(z)] 



(5) 



is the symmetrically ordered characteristic function of the 
state p. A state is called Gaussian if its characteristic func- 
tion is a Gaussian function in phase space A Gaus- 
sian map is a completely positive map that maps all unknown 
Gaussian states onto Gaussian states and a Gaussian unitary 
a unitary generated by a quadratic polynomial in the canoni- 
cal coordinates, reflected by a symplectic transformation from 
Sp(2n, R) on the level of canonical coordinates. 

In the construction of the Gaussian unitary representa- 
tions ([TJ of $, the vector v plays a marginal role since it can 
be eliminated via a unitary rotation acting on the output state, 
see, e.g., Ref. B1411 . In contrast, the matrices in Eq. OJ are 
of fundamental importance — in particular, we shall see that 



the value of £p*L and of our estimation of depend upon 
the ranks of Y, £, and Y — iE. It is thus worth anticipating 



*(*) 



some relevant facts that concern these matrices. First of all we 
notice that the inequality @ implies the following relations 

ker[E] n ker[F - tS] C ker[y] C ker[F - jE] , (6) 
ker[F] C ker[E] , (7) 

where throughout the paper, given a generic (possibly real) 
d x d matrix M, we denote its kernel (null subspace) with 
ker[M] := {w e C d : Mw = 0} [20]. The first inclusion 
in Eq. (0 follows from the definition, the remaining one and 
the inclusion of Eq. (0 are derived from the observation that 
w^Yw = => w^{iZ)w = => w^(Y - iE)w = => 
(Y — iE)w = Ew = 0. Putting these identities together 
we also find that 

ker[F] = ker[E] n ker[F - iE] . (8) 

Other useful properties are the identities 

2 rank[F - iE] = rank[F] + rank[F - EF el E T ] , (9) 

and the inequalities 

rank[F] ^ rank[E] > rank[F] - rank[F - EF el E T ] > , 

(10) 

where rank[M] stands for the rank of the matrix M (i.e. the 
dimension over the complex field of the complement to C d of 
the matrix kerfM]), and Y el is the Moore-Penrose (MP) in- 
verse of Y 112111 . The explicit proof of these relations is rather 
technical and thus we postpone it to Appendix lAl Here we 
rather point out that the first inequality of Eq. ( fTOb is a conse- 
quence of the fact that ker[Y] is included in ker[E], while the 
last inequality is an immediate consequence of the fact that 
EF el E T is positive semi-definite. 

In Ref. fil, an upper bound for was set by showing 
that one can construct a Stinespring dilation of $ that involves 
I = 2n — r'/2 environmental modes with 



:= rank[F] - rank[F - Er el E T ] 



(ID 



In what follows we will strengthen this result by showing that 
the minimum number of modes necessary to build a Gaussian 
Stinespring unitary dilation for $ is given by 



^ = rank[y] - r'/2 = rank[y 



iE] 



(12) 



where we used Eq. (O when formulating the second identity. 
Since Y is a 2n x 2n matrix, we have 2n — k ^ 0, and so the 
optimal bound we prove here leads to a significant improve- 
ment compared to the results of Ref. (3. In particular, for 
those BGCs which represent unitary transformations of the n 
input modes (i.e., Y — and X e Sp(2n, R) symplectic lU]) 
the optimal bound ( fT2b yields £ p*r e = — no environment is 
required to construct the dilation — while Ref. flill had this 
value equal to 2n. To prove Eq. dT2b we shall first show that 
the quantity rank[Y — i£] provides a lower bound for £p2 e (see 
SecHHI) and then construct an explicit Stinespring dilation (Q]) 
for $ which attains such bound (see Sec. lIVI ). 
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III. 



LOWER BOUND ON 4*re VIA GENERALIZED 
CJ-STATES OF BGCS 



In this section we review the notion of generalized Choi- 
Jamiolkowski (CJ) state for a multi-mode BGC (see also 
Ref. lfl8ll and compare Refs. lfl9ll ). and use it to show that the 
term on the rhs. of Eq. ( fT2] i provides a lower bound for £ptre- 
Consider a state vector \^^)a,b providing a purification of a 

quantum state A = J^jLi A j(\j){j\)A of the system labeled 
A which has full rank (e.g., a Gibbs state of n modes). That is 
to say, 



l*A> 



A.B 



E 

j=0 



= (A^®t)J2\j)A®\j}B, 
3=0 

with A indicating the input space of the channel $, B being an 
ancillary system isomorphic to A, and {\j) : j = 0, • • • , oo} 
denoting an orthonormal complete basis. A generalized CJ 
state of the channel $ is now obtained as 

fct,fl(S) = (*®Z)(|*A><*AlUBi ( 13 > 

with I being the identity map. The state pa,b(&) provides a 
complete representation of the channel via the inversion for- 
mula, 

$(p) = Tr B [(lU ® A B 1/2 plA B 1/2 )p A , B m , (14) 

where ps and A g are copies of the states p and A on B, re- 
spectively, while pig is its transpose with respect to the or- 
thonormal basis introduced above. We will suppress an index 
labeling both the chosen basis and the reference state. 

For finite-dimensional system pa,b($) provides a standard 
CJ state representation when A is taken to be the maximally 
mixed state (compare Refs. lfl9ll ). In the infinite-dimensional 
case such limit in general is well defined only in the context 
of positive forms, see Ref. 1181 . However, Eq. ( fT~4T > shows 
that we do not need to approach such a limit in order to build 
a proper representation of the channel: It is defined for any 
state diagonal in the distinguished basis of full rank. Further- 
more, it is easy to verify that it is always possible to work with 
CJ states pa,b(®) which are Gaussian: To do so simply take 
(l^s) {9^}) a.b to be Gaussian and use the fact that the Gaus- 
sian map $ ® I maps Gaussian states into Gaussian states. 
In the following we choose to take such Gaussian reference 
states. In particular, we will assume ^) ("^ t\) a,b to be a 
Gaussian purification of a multi-mode Gibbs (thermal) state 
of quantum mechanical oscillators. 

An important observation concerning the generalized CJ 
representation is that, given a Stinespring representation of $ 
involving an environmental system E, one can construct a pu- 
rification of pa,b{$) that uses E as ancillary system. Indeed, 
assuming that Ua.e and (|0)(0|)b give rise to a Stinespring 
representation for we have that the pure state with state 
vector 



\x)a.b,e = U a ,e\*a) 



A.B 



® |0)j 



is a purification of pa.b{^)- Furthermore, if pa,b(&) is 
Gaussian and E represents a collection of £ environmental 
bosonic modes with \0}e being a Gaussian state vector and 
Ua,e being a Gaussian unitary, it follows that also \x)a,b.e 
will define a Gaussian purification of pa,b(&)- Putting these 
facts together it follows that a lower bound for the minimal 
number ^p*e of environmental modes that are needed to build 
a Gaussian Stinespring representation of <I> is provided by the 
minimal number q^} of Gaussian ancillary modes that are re- 
quired to purify a generalized Gaussian CJ state pa,b(&) of 
$, i.e., we have that 



/(*) > 

'-pure f 



„(*) 



(16) 



To compute gffi we first make a specific choice for |^ i)a,b- 
In particular, since A is composed by n bosonic modes, we 
can take I^^a.b to be a product of n identical two-mode 
state vectors of the form 



\*x)A,B = 6d\Tp)A i ,B i 



(17) 



where reflects a purification of a Gibbs state of the 

i-th mode A; of A built by coupling it with the correspond- 
ing ancillary system Be This is nothing but what is usually 
referred to as a two-mode squeezed state lfl5ll . The resulting 
state is of course Gaussian and it is fully characterized by its 
covariance matrix. To express it in a compact form note that 
the kernel of the natural symplectic form for the 2n modes of 
A, B is given by 



o~ A,B := 



T2n 








C2n 



(18) 



where the upper-left and lower-right block matrices represent 
the symplectic forms of the n modes of A and B, respectively, 
defined as in Eq. (0J. With this choice the covariance matrix 
7 of (\^ ^)a,b(^ \\) is given by the following GZ(4n, R) ma- 
trix, 



7 



a 


5 ' 


8 T 





(19) 



where a, j3 € Gl(2n, R) are the covariance matrices of the A 
and B modes, respectively, with S, 5 T G Gl(2n, R) being the 
cross-correlation terms. Explicitly, they are given by 



a 














with 8 > 1 and 



f{6) := -(0 2 - l) 1 / 2 . 






f(d)t n 


J(0)t n 






(20) 



(15) 



The parameter 8 determines the temperature of the Gibbs 
states we used to build the vector \*&^)a,b< or equivalently, 
the two-mode squeezing parameter of the purification. In par- 
ticular, the case 8=1 corresponds to the limit in which all 
the modes of A and B are prepared into the vacuum state: In 
this case the state A no longer has maximum support and thus 
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does not provide a proper starting point to build a CJ state. For 
9 — > oo, in contrast, the state \^^)a,b approaches a purifica- 
tion of a maximally mixed state for the modes (for details see 
Ref. (HI]). Equivalently, it corresponds to the limit of large 
squeezing in the two-mode squeezed state of the purification. 
Notice also that by construction, for all values of 9 ^ 1, 7 
satisfies the condition 7 ^ iua.b, as it indeed represents a 
physical pure state. 

The generalized CJ state pa.b(®) for a Gaussian channel 
characterized by matrices Y and X as in Eq. (0 is now com- 
puted as in Eq. (FX 3t . The resulting state is still Gaussian and 
has the covariance matrix 7' E Gl(An, R) given by 



7 = 
where 



X T aX + Y 


X T S~ 




ex T x + Y 


f(8)X T a x ' 


6 T X 






f(6)a x X 


Bt 2n 






In " 


In 






(21) 



In general it will be a mixed state and we are interested in the 
minimum number q^ n of ancillary modes q that is needed to 
construct a Gaussian purification of it. As discussed in Ap- 
pendix|Bj this is given by the quantity 



,(*) 



rank[7' - ia A ,B] 
2n — dimker[7' - 



- 2n 
ioA.s\ , 



(note that in this case 7', <ta,b & GZ(4n x in, R)). In what 
follows we will compute this quantity, showing that it coin- 
cides with the right hand side of Eq. (fT2l . To do so, we first 
notice that the dimension of the kernel of 7' — iua,b can be 
expressed as 



dim ker [7' —i<JA,B\ — dimker 



dX T X + Y -ia 


f(9)X T 


f(6)X 


0t-2n + icr 



(23) 

where the second identity was obtained by rotating 7' — wa.b 
with the transformation 



T := 











Ox 



(24) 



As for any positive semi-definite matrix M, the kernel in 
Eq. d23l can be computed as the set of vectors w G C d 
which satisfy the condition w^Mw = lf20ll . Writing w — 
(wa, wb), we arrive at the condition 

9 (w* A X T Xwa — w*aX t wb — wgXwA + w* b wb) 
+ w* A (Y — ia) wa + WgiawB 
+ 0(1/9) (w A X T w B + w* B Xw A ) = 0, (25) 

where in the first and second line we have collected all terms 
which are linear and constant in 9, respectively. For 9 > 1 
sufficiently large this requires the following conditions, 

w* A X T Xwa — w* a X t wb — w* b Xwa 

+w* B w B = , (26) 
w* A (Y — ia) wa + w B io"WB = . (27) 



The first equation means Xwa = vjb, whereas the second 
reads w* A (Y — ia) w A + w* A iX T ctXwa = 0, that is 



w* A [Y -iE}w A = 0. 



(28) 



There is one-to-one correspondence between solutions wa of 
Eq. d28b and w = (wa, Xwa) of Eq. d25l ), hence 

dimker[7' — icta.b] — dimker[Y — i£]. 

Replacing this into Eq. d22b we finally get 

q^l = 2n- dimker[F - iS] = rank[y - iE] , (29) 

where in the last identity we used the fact that Y — iH is a 
2n x 2n matrix. 

IV. OPTIMAL BOUND AND EXPLICIT CONSTRUCTION 

In this section we explicitly construct a Gaussian unitary 
dilation with q^) n = rank[Y — iH] environmental modes. In 
this way, we demonstrate that the lower bound derived in the 
previous section is tight, concluding the derivation of Eq. ( fT2b . 
To do so, let us assume that the number of modes which define 
the state pe in Eq. (Q]) are q^ n . Without loss of generality, we 

write the kernel of the form corresponding to the commutation 

(<&) 



(22) relations of our n + q^^ modes in block structure 



&2n 



02n 











(30) 



where oi n and a E m are 2n x 2n and 2q^ n 



x 2q 



ma- 



trices associated with the system and environment, respec- 
tively. While (T2n is defined as in Eq. ([3), for a E (s) we do not 

make any assumption at this point, leaving open the possibil- 
ity of defining it later on. Accordingly, the Gaussian unitary 
IJ A e of Eq. ([T|i will be determined by a symplectic matrix 



S e Sp(2(n + q^l), R) of block form 



S := 



si 


S2 


S3 


S4 _ 



(31) 



satisfying the condition ScrS T = a. In the above expressions, 
s\ and S4 are 2n x 2n and 2q^ n x 2q^ n real square matrices, 

while S2 and sj are 2n x 2q^ n real rectangular matrices. As 
noticed in Ref. 03], the possibility of realizing the unitary 
dilation (Q~|i can now be proven by simply taking 



si = X 1 



(32) 



and finding S2 and a g^^-mode covariance matrix [1] je ^ 



10 



2qW 



satisfying the conditions 



S2lES 2 = Y, 



(33) 



with -fE being the covariance matrix of the Gaussian state p B 
ofEq.fl}. 
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First, let us consider the case in which q^y is an even num- 
ber. To identify valid s 2 and 7^ which solve Eq. ( l33l >, it is 
useful to transform Y and E as in Eq. ( |A7t and dA8t of Ap- 
pendix [A] (take A = T, B = E, m = 2n, a = fc, and 
b = r = rank[E]). Actually, applying an extra orthogonal 
matrix, Y' is still like in (1A71 >, while E' can be written as 



E' := C*EC T = 






M 









































•r)/2 



}r/2 
}(,<*> 
}r/2 (34) 

}hSl-r)/2 

}2n ■ 



(*) 

'mm ' 



where C G Gi(2n,H) and fj, = diag(/xi, • • • ,fJ, r /2) is the 
r/2 x r/2 diagonal matrix formed by the strictly positive 
eigenvalues of |E'| (satisfying i r / 2 /U as in Appendix lAb. 
Introducing then s' 2 := C s 2 the conditions of Eqs. fl33l can 
be equivalently written as 



s'^e{s' 2 Y =Y' . (35) 



The explicit expressions for corresponding 7^ and s 2 are ob- 
tained in the following way. We take the environmental sym- 
plectic form to be 



Ok © tTfc_ 



(36) 



where we have set A; :— rank[y]. A unitary dilation with 
#miri = k — r'/2 environmental modes in a pure state is ob- 
tained by choosing the 2n x 2q^ n rectangular matrix s' 2 as 

(37) 



' K- 1 


A ' 









with K being the fc x fc symmetric matrix defined by 



K :-- 














l(fc-r)/2 





(U-V2 








l(fc-r)/2 



(38) 



and A being a rectangular matrix k x (A: — r') of the form 



A := 













}r'/2 















}(r- 


r')/2 









l(fc-r)/2 


}fc/2 


- r/2 












}r'/2 















}(r- 


r')/2 





l(fc-r)/2 






}fc/2 


- r/2. 



(39) 



By direct substitution one can easily verify that the first con- 
dition of Eq. (l35T l is indeed satisfied. Vice versa, expressing 
the (2k — r') x (2k — r') covariance matrix of /5g as 



IE 



a 


8 ' 


8 T 





(40) 



the second condition of Eq. (135 D yields the following equation 

a + A5 T + 5A T + A/3A T = K 1 . (41) 



A solution can be easily derived by taking the k x k block a 

as 



























5 11 

4 J1 (fc-r)/2 



(42) 



while /3 and S are, respectively, (fc — r') x (fc — r') and k x 
(/c — r') real matrices defined as follows: 



/?:= 



Mo 1 











fl(fc-r)/2 





Mo 1 








|l(fe-r)/2 . 



, (43) 



S := 


















/(Mo" 1 ) 












3 11 

-4-U(fc-r)/2 












/(Mo" 1 ) 













-|l(fc-r)/2 







(44) 



with /i D is the (r — r') /2 x (r — r')/2 diagonal matrix formed 
by the elements of /i which are strictly smaller than 1, and 
with f(9) defined as in Sec. [HI] Notice that the parameter r' 
(defined above) corresponds also to the number of eigenvalues 
having modulus 1 of the matrix E', i.e., 



2n — rankfls 



E'(E') 



(45) 



as can be easily shown by using Eq. JA18t with A = Y and 
B = E. With the choice we made on the commutation matrix 
(T , the matrix a is a fc x fc covariance matrix for a set of 

2 ?min 

independent fc/2 bosonic modes, the matrix f3 is a (fc — r') x 
(k — r') covariance matrix for a set of independent (fc — r')/2 
modes, and the matrices 8 and S T represent cross-correlation 
terms among such sets. For all diagonal matrices /i compatible 
with the constraint 



II 



r/2 ^ A* ) 



(46) 



the solution je satisfies also the uncertainty relation je ^ 
ia E . Furthermore, since it has 

Det[ 7£ ] = 1, (47) 

this is also a minimal uncertainty state, i.e., a pure Gaussian 
state of q^ n modes By a close inspection of the covari- 
ance matrix 7^ derived above, one realizes that it is composed 
of three independent pieces. The first one describes a collec- 
tion of 7-'/2 vacuum states. The second one, in turn, describes 
(r — r')/2 thermal states characterized by the matrices /Lt" 1 
which have been purified by adding further (r — r')/2 modes. 
The third one, finally, reflects a collection of fc — r modes pre- 
pared in a pure state formed by fc/2 — r/2 independent pairs 
of modes which are entangled. Let us point out again that 
this covariance matrix is indeed formed by q^ n modes. The 

whole derivation can be trivially extended for q[^ n odd, by 
adding to the previous covariance matrix a single mode in the 
vacuum state. 
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V. DILATIONS WITH MIXED ENVIRONMENTS 



real matrices: 



In Ref. II 1-411 it was shown that for arbitrary (not necessar- 
ily Stinespring) dilations one can consider an environment of 
only I = 2n — r/2 modes — observe that r is larger than the 
quantity r' introduced in Sec. HI] because of Eq. (TTOb . Here, 
we will strengthen this bound by showing that it is possible to 
construct a unitary dilation using just 



- k 



r/2 = rank[Y] -rank[E]/2 



(48) 



environmental modes which are prepared in a Gaussian, but 
not necessarily pure, state. Note that the term on the rhs. is 
nonnegative due to the first of the inequalities in Eq. ( fTUl ), and 
that it is explicitly smaller than the one provided in Ref. 11411 
due to the fact that Y is a 2n x 2n matrix. It is worth stressing 
however that differently from the pure dilation case, we are 
not able to determine whether Eq. d48l l is indeed the optimal 
bound (we believe it is). 

For the sake of simplicity, again we will treat explicitly only 
the case of k even (the analysis however can be easily ex- 
tended to the odd case). Because of the structure of A given in 
Eq. d39l ), the (k — r) environmental modes prepared in a pure 
state (see the end of Sec. |IV| ) enter explicitly in the identity in 
Eq. (HTb : consequently, if we wish to satisfy such relation, we 
cannot remove any of these modes without changing A. Vice 
versa we can drop some of the auxiliary modes which were 
introduced only for purifying the environmental state. Since 
they are (r — r')/2, we can reduce the number of modes from 



£ W to 



#(*) 



= £$>-(r-r')/2 = k-r/2. 
To see this explicitly, take 



re 



CTfc © <Xfc_ 



(49) 



(50) 



The matrix s' 2 can be still expressed as above but with A being 
a rectangular matrix k x (fc — r) of the form 













A := 


l(fc-r)/2 


(51) 










*(*— r)/2 





Similarly, f3 and 5 entering in the definition of je become, 
respectively, the following (fc — r) x (fc — r) and k x (fc — r) 



/?:= 



5 11 








|l(fc-r)/2 . 



(52) 



and 












3 11 














(53) 



This covariance matrix now consists of two independent parts: 
the first one describes a collection of r/2 thermal states de- 
scribed by the matrices /i -1 . The second one reflects a col- 
lection of k — r modes prepared in a pure state formed by 
k/2 — r/2 independent couples of modes which are entangled. 
VI. CONCLUSIONS 

We have analytically computed the minimum number of en- 
vironmental modes necessary for a Gaussian unitary dilation 
of a generic multi-mode bosonic Gaussian channel. Moreover, 
we have also explicitly demonstrated how to construct such 
a Gaussian dilation in terms of the covariance matrix of the 
noisy environment and the symplectic transformation associ- 
ated to the unitary system-environment interaction. These re- 
sults may allow one to introduce a classification of the bosonic 
Gaussian channels in terms of the corresponding noise in- 
duced by these maps, which is somehow related to the mini- 
mum number of environemntal modes to represent such chan- 
nels. Moreover, constructing a dilation with a minimal num- 
ber of auxiliary modes may be useful to minimize the size of 
the corresponding complementary channel and then to sim- 
plify the degradability analysis, which is extremely useful in 
the calculation of the quantum capacity of these continuous- 
variable quantum maps. 
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Appendix A: An important identity 

In this Appendix, we prove the important identity (O and 
the inequality dlOl l. by the following more general lemma. 

Lemma 1 Let A, B <E Gl{m, R) be mxm real matrices with 
B being skew-symmetric, which satisfy the inequality 

A^iB . (Al) 

Then given A el the MP inverse ^2~/l/ of A, the following iden- 
tity holds 

2 rank[A - iB] = rank[A] + rank[A - BA el B T 
Furthermore the following inequality applies 



■(A2) 



rank[B] > rank[A] - rank[,4 - BA B B ] . (A3) 

Proof: Let us start by reviewing some general properties of 
A and B. Because of Eq. dAU the matrix A must be positive 
semi-definite, and its support must contain the support of B. 
Consequently indicating with a = mnk[A] and 6 = rank[B] 
the ranks of the two matrices, we must have a ^ 6 with b even. 



Furthermore, defining IT S Gl(m, R) to be the projector on 
the support of A, it will commute with A and B and hence 
satisfy the following identity 



nA = An = A, iib 

Consider then the invertible matrix 

A:= A+{l m -U) 
The MP inverse B22I1 of A is defined by 
A el IL-l- 1 !! . 



B n = B . (A4) 



(A5) 



(A6) 



To prove the validity of Eq. ( IA2b we note that it is possible 
to identify a congruent transformation A i-> A' = CAC T , 
B^tB' = CBC T , with C G Gl(m, R) invertible such that, 



A' = 








}a 








}m — a 



(A7) 



and 



B' = 













}b/2 


-A« 








}b/2 










}a-b 








} m — a 



(A8) 



with fi = diag(/ii, /i2, ■ ■ ■ ,A t &/2) being the 6/2 x 6/2 diago- 
nal matrix formed by the strictly positive eigenvalues of \B'\ 
(by construction they satisfy 1 ^ fij ^ 0). The matrix C can 
be explicitly constructed as follows. First we identify the or- 
thogonal matrix O 6 Gl(m, R) which diagonalizes A and IT 
puts them in the following block forms: 



OAO T = 



ono T 



A" 





}a 








}m — a 


" la 





}a 








} ra — a 



(A9) 



(A10) 



with A" G Gl(a, R) being a a x a positive definite diag- 
onal matrix. Then we construct the invertible matrix K G 
Gl(m, R) defined as 



K = 








}a 





1 m a 


} jri — a 



(All) 



(notice that the matrix A"~ 1/2 G Gl(a, R) is well defined 
since A" G Gl(a, R) is invertible). Finally, we take O' G 
G/(a,R) to be an orthogonal a x a matrix and define C as 
follows 



C = 



' O' 








1 m — a 



KO = 



' 0'A"~ 1/2 








1 m — a 



O .(A12) 



By construction we have that for all the choices of O' the 
resulting matrix is invertible and Eq. ( IA7b is satisfied. Vice 
versa, Eq. ( IA8b can be satisfied by noticing that, since the sup- 
port of B is included into the support of A, we must have 



KO B T K T = 



' B" 





}a 








}m — a 



(A13) 
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with B" £ Gl(a, It) skew-symmetric and having the same 
rank as B. By using a theorem from linear algebra one can 
then find an orthhogonal O' £ G7(a, It) such that 



/T 























(A14) 



with p, being a positive diagonal matrix of dimension equal to 
the rank of B" (the elements ±ipj are its not null eigenval- 
ues). Using such an O' in order to build C as in Eq. ( IA12I ) we 
can then satisfy Eq. ( IA8b . 

Now we notice that, since any congruent transformation 
preserves the rank of a matrix, the following identity holds: 

rank[A - iB] = rank[C(A - iB)C T ] 







—ifi 







rank 

























where counts the number of eigenvalues of the matrix 

/i which are equal to 1. The last identity follows from counting 
the non-zero eigenvalue of the matrix on the left-hand-side of 
the second line. This can be easily done by observing that its 
spectrum contains m — a explicit zeros (these are the terms 
in the zero block diagonal term), a — b ones (these are the 
ones on the diagonals of the first block) and 1 ± p,j with fij £ 
[1, 0] being the eigenvalues of p. Consequently, the non-zero 
eigenvalues are obtained by subtracting from k (rank of the 
first block) the number #i(/x) of eigenvalues of p which are 
equal to 1. To compute the latter quantity we note that 



B'B' T = 



\ 










}b 










}b 










}a-b 








}m — a 



(A15) 



which yields 

rank[l m - B'B' T ] =m-2# 1 (/x) . (A16) 

Using the fact that C is invertible, one has 

rank[l m - B'B ,T ] = rank[l 7 

= rank[C* _1 C~ T 

Since O' and O are orthogonal, we notice that C~ 1 C~ T is 
composed of two terms that span orthogonal supports. Specif- 
ically we can rewrite it as 



CBC T CB T C T ] 
BC T CB T \ . 



C- X C- T = O t KO = T 



A + O 7 












m — a 



o 



' A" 








1 m — a 


A + 


(l m - 



o 



n) = a, 



where Eqs. (|A9t and (1A1 Q have been used. Similarly, 
BC T CB T is only supported on the support of A. Indeed, 
we have 

bc t cb t = (n b n) [c- 1 ^]- 1 (n b t n) 
= (n b n) A- 1 (n b t n) = (n b n 2 ) A- 1 (n 2 b t n) 
= (n b n)(n A' 1 n)(n b t n) 
= (n b u)A el (n b t n) = b a &1 b t . (An) 



Using these identities, we can then rewrite Eq. ( IA17I ) as 

rank[l m - B' B' T \ = rank[A - BA G1 B T ] 

= rank[l m - n] + rank[A - BA el B T ] 

= m-a + rank[A - BA el B T ] . (A18) 

Thanks to Eq. (IA16b the above identity finally yields 

rank[yl] - rank[A - BA el B Tl 



(A19) 



which gives Eq. (1A2I) when inserted into Eq. dA16l ). The in- 
equality (IA3t can finally be proven by noticing that because 
of the invertibility of C, one has rank[_B] = rank[_B'] = b 
which, by construction, is larger than 2#i (/i). The result then 
follows simply by applying Eq. ( IA19I ). ■ 



Appendix B: Gaussian purifications 

Here, we emphasize the minimal number of ancillary 
modes which are necessary to construct a Gaussian purifica- 
tion (i.e., a purification which is joint pure Gaussian state of 
the system and of the ancillary modes) of a generic multi- 
mode Gaussian state p. Of course, the Gaussian requirement 
on the purification is fundamental for our purposes: Since any 
number of modes can always be embedded in a single one, by 
dropping it the minimal number of modes is always smaller 
than or equal to one. 



1. Minimal Gaussian purifications of Gaussian mixed states 

Let 7 £ Gl(2n,TR) the covariance matrix of a Gaussian 
state p of a system A formed by n bosonic modes. We know 
that it must satisfy the following inequality 



7 ^ i°2n 



(Bl) 



with <j2n £ Gl(2n,WL) being the skew-symmetric matrix 
in Eq. © representing the symplectic form of the modes. 
Thanks to the Williamson's theorem H23I1 we know that there 
exists a symplectic transformation S £ Gl(2n, R) which al- 
lows us to diagonalize 7 in the following form 



7 1 — ^ S'yS 1 



D 








D 



(B2) 



with D £ Gl(n,~Rl) being the diagonal matrix formed by 
the symplectic eigenvalues Dj of 7 which satisfy the con- 
dition Dj ^ 1 as follows from Eq. ( 1BU . The values {Dj} 
are the symplectic eigenvalues of 7 (H Hal , so the positive 
square roots of the eigenvalues of the matrix —<J2nl&2nl £ 
Gl(2n,WC). The transformation 7 i-> S^fS T corresponds to 
appling a Gaussian unitary to the state which transforms it 
into a product state of the n modes, in fact a product of Gibbs 
states of unit harmonic oscillators. Hence, it does not restrict 
generality to assume that 7 is of the form of the rhs. of Eq. 
( IB2I ) in the first place. 

Now, let r be the covariance matrix of the minimum pu- 
rification of 7, viewed as being defined on a bi-partite system 
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labeled A — the original system — and B. Since the spec- 
trum of the reduced state with respect to B is identical to the 
spectrum of the reduced state of A, also the symplectic spectra 
of the two reductions are the same. Hence, it does not restrict 
generality to take T to be of the form 



D 





c 





D 


c 


T 


D 










D 



(B3) 



with suitable C € Gl(2n,M) such that the symplectic spec- 
trum of r consists of 1 only, with respect to the symplectic 
form in the convention as in Eq. (fT8T l. Now, by taking 



C = 






'/ 








(B4) 



with r\ = diag(/(.Di), • ■ • ,f(D n )), one clearly arrives at 
the covariance matrix of a valid purification. This purifica- 
tion essentially involves as many modes as there are symplec- 
tic eigenvalues different from 1 — those modes associated 



with unit symplectic eigenvalues correspond to pure Gaus- 
sian states already. Denoting the number of unit values in D 
by this purification hence involves n — #i(D) many 

modes. Invoking the definition of the symplectic spectrum, 
one finds that 



rank[ 7 - a 2n l cr 2 J/ 2 



(B5) 



It is also easy to see, however, that no purification can involve 
fewer modes than that. Consequently we have 



(B6) 



The covariance matrix of the reduced Gaussian state of the 
purification with respect to B is necessarily given by the 
rhs. of Eq. (IB2I) . up to local symplectic transformations S € 
Sp(2n, R). Hence, any Gaussian purification must involve at 
least involve n — #\(D) modes, as so many symplectic eigen- 
values are different from 1. Needless to say, if one gives up 
the property of requiring a Gaussian purification, one can al- 
ways embed the purification in a single mode, if the state is 
mixed, while no additional mode being required if the state is 
already pure. 



